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Summary 

A short self contained computer program for calculating sky-wave interference 
at mf. has been written. The program employs a simple empirical law for calculating 
field-strength and makes a correction for the transmitting aerial vertical radiation pattern, 
the take-off angle being derived from the great circle distance and the number of 'hops'. 
Protected field-strengths are calculated for each wanted service area from the interference 
levels from individual co-channel and adjacent-channel sources, from all co-channel 
transmitters taken together, from all adjacent- channel transmitters taken together and 
from all co- and adjacent-channel transmitters taken together. 
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1. Introduction 

This report describes a short program to compute m.f. 
sky-wave interference which has been designed as an aid in 
the preparation of international plans for m.f. broadcasting. 
Because it uses recognised prediction methods and is very 
economical to run, it could provide a ready means for 
testing a number of possible plans or be used to assist in the 
assessment of ad hoc modifications to these plans. For the 
calculations the transmitter information is read in from 
cards, permitting maximum flexibility in setting up new 
plans. The cards hold details of all the known transmitters 
to be considered from which any selection may be made to 
run the program. If required, the frequency of any station 
may be changed and various options are available which are 
explained in Section 3. 

The purpose of the program is to assess the effect of 
co-channel and adjacent-channel interference by sky-wave 
propagation in service areas centered on any number of 
specified transmitters. (These are called the wanted trans- 
mitters). The program selects those transmitters which are 
co-channel and adjacent-channel with each wanted trans- 
mitter in turn and computes protected field-strengths based 
on the sky-wave interference from these interfering sources, 
separately and in combination. The combined effect of 
interfering sources is derived by taking a weighted sum of 
the interfering signal powers, the weighting being in accor- 
dance with protection ratios which are determined by the 
frequency spacing between the wanted and interfering 
carriers. 

A distinction is made between interfering stations 
which are on precisely the same frequency as the wanted 
signal as to whether they will be carrying the same or 
different programme material. This affects the choice of 
the protection ratio. 



2. The input cards 

The transmitter input cards contain the following 
information:— 



(i) Frequency (a whole number of kHz) 

(ii) Station number (numbered in frequency sets) 

(iii) Station name 

(iv) Country 

(v) Latitude and longitude 

(vi) Transmitter power (in dB rel. 1 kW) 



(vii) Programme code (this enables the program to dis- 
tinguish between stations carrying the same or dif- 
ferent programme material). 

A set of change cards can be added to the transmitter 
cards to indicate frequency changes to selected stations. 
These additional cards also permit options on the selection 
of wanted and interfering transmitters by the program. 
The options are described in the next section. 



3. Options 

There are six options for running the program; five of 
these require that the option variable is set to the appro- 
priate value on the change cards. If no change cards are 
included in the input the program will run in its basic form 
in which each station from the input set is treated as wanted 
in turn, co-channel and adjacent-channel stations are then 
selected from the remainder and treated as interfering 
sou rces. 

3.1. 'Cards' 

When the option variable is set to 'CARDS' the total 
set of transmitters considered is restricted to those appear- 
ing in the change cards. Otherwise the program runs as 
above (frequency changes may, of course, be to the same 
frequency if no alteration is required). 

3.2. 'Single' 

In this case only those stations listed by the change 
cards are considered as wanted; all stations which are input 
on existing sets or change cards are used as potential 
sources of interference. 

3.3. 'Europe', 'Africa' or 'Asia' 

If the option variable is set to either of these con- 
tinents wanted transmitters are only taken from the con- 
tinent indicated. 



4. Vertical radiation pattern (VRP) 

Three types of vertical mast-radiator aerials are provided 
for. These are 'short', 'A/4' and 'A/2' respectively. The 
v.r.p. corrections applied for these types follow the recog- 
nised EBU method 1 and are shown in Fig. 1. The correc- 
tion to be used is chosen in accordance with the specified 
transmitter power. 
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elevation angle, degrees 

Fig. 1 - Vertical radiation patterns for transmitting aerials 

normalised to c.m.f. = 300 Volts 

(fl) Short aerial (b) A/4 aerial (c) A/2 aerial 

Type 1, 'short', transmitter power <1 kW 

Type 2, 'A/4', transmitter power >1 kW, <50 kW 

Type 3, 'A/2', transmitter power > 50 kW 

This procedure has been adopted as the actual aerial types 
are not known in the majority of cases. In practice aerial 
heights tend to increase with transmitter power roughly in 
this way. These curves are normalised to a c.m.f. of 300 
Volts which is nominally that of a 1 kW transmitter and a 
short aerial (100% efficient) in the horizontal direction. 
The appropriate correction is applied by assuming that the 
reflection point in the ionosphere is 100 km above the 
surface of the earth and calculating the take-off angle also 
taking into account the minimum number of 'hops' which 
are involved. The minimum number of 'hops' for each 
path is determined from the great circle distance. The 
radiation patterns are assumed to be identical on all 
bearings. 

The v.r.p.'s take no account of the reduction of field 
strength at very small angles of elevation due to poor ground 
conductivity close to the aerial. This fact may introduce 
error in a small percentage of the field strength calculations 
but will not be serious for planning estimates. 



5. The median field-strength 

The median field-strength (dB rel. 1/uV/m) predicted for 
each transmitter in the wanted service area is calculated 



from a simple experimental* formula which is a function of 
the great circle distance only. This formula is known to be 
reasonably accurate for distances greater than 3000 km. 
Referring to Fig. 2 it is seen that it also corresponds closely 
to the CCIR curve for distances less than 2500 km taken 
at a frequency of 1000 kHz. The frequency dependence 
of the CCIR estimate is weak and it appears that the pro- 
posed curve is very adequate at all distances for planning 
purposes. 

A flow diagram for the field-strength sub-routine is 
given in Appendix III. 



6. Interference weighting 

The program specification requires that the protected 
field-strength in the service area should be calculated in view 
of each interfering transmitter taken individually and in 
view of all co-channel transmitters taken together and all 
adjacent-channel transmitters taken together. The com- 
bined effect of many interference sources is achieved by 
weighted signal power addition, the weighting curve of 
Fig. 3 being used for the purpose.t As specified, fre- 
quency spacings of up to 4 kHz cover co-channel sources; 
spacings from 5-10 kHz cover adjacent-channel sources. A 
separate curve based on frequency spacings can be em- 
ployed when it is necessary to distinguish between trans- 
mitters carrying the same or different programme material. 
In the present case the 'same programme' situation only 



* Contained in an interim document of the EBU which is not yet 
published. 

t The protection ratios are being reconsidered and may be modified 
in the future. 
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Fig. 2 - Median field-strength with great circle distance for 

cm. f. = 300 Volts, X = 300 m 
(1) Full line = 80-2 - 10 |og 10 D - 0-00777 X"" 0-26 D dB, 
CCIR curve 
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empirical law used in the program 
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Weighting curve (protection ratios) with frequency 
difference 



applies when the interfering source has zero frequency 
spacing. The weighting curve is included in the program as 
a data statement. No interpolation is required as only 
integral values of the argument are allowed. 



7. Output 

An example of line printer output isgiven in Appendix II. 
The output data includes, amongst other information: 

(1) the name of the transmitter serving the wanted area; 
its frequency and code number: 

(2) a list of all interfering transmitters which are either 
co-channel or adjacent-channel, their existing and 
newly assigned frequencies, the protected field- 
strength required in view of each one taken alone and 
the appropriate protection ratios: 



(3) the protected field-strength required in view of all the 
co-channel interfering signals taken together: 

(4) the protected field-strength required in view of all the 
adjacent-channel interfering signals taken together: 

(5) the protected field-strength required in view of all the 
co-channel and adjacent-channel interfering signals 
taken together. 

At the end of each run the number of equal protected fields 
for multiple co-channel, adjacent-channel and the combined 
cases are totalled in dB steps. 

8. Discussion 

The computer program described in the foregoing is 
intended as a tool for m.f. service planning on an inter- 
national basis. More sophisticated programs are being 
considered for the future and these would be expected to 
contain extensive data banks of, for example, ground 
conductivity, magnetic dip angles, etc. It should be 
remembered, however, that computer data banks can be 
costly to store and handle in a program and the probability 
that the cost and complexity will be justified in terms of 
improvements in prediction will have to be carefully 
assessed. The short self-contained program described can 
easily be adapted for use on most computers and, therefore, 
it might be arranged to be available at international meetings 
as an on-the-spot aid to negotiations. The only sub- 
routines required are those for standard mathematical 
functions. 
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Appendix I 

INPUT DATA (Cards) 
A- EXISTING STATIONS 



Variable 


Forma 


Name 




NTXS 


14 


KHZCOD 


16 


NAMCUN 


4A6 


ALAT 


14 




13 




A1 


ALONG 


14 




13 




A1 


POWDB 


13 


ISYNC 


13 



Total number of transmitters in existing set 
Frequency and station number 
Name and country 

Latitude in degrees and minutes (N or S) 

Longitude in degrees and minutes (E or W) 

Transmitter power in dBs 
Programme code 



B- CHANGE CARDS 



NFRE 
NUMNEW 

IWAREA 



14 New frequency 

12 No. of transmitters to be allocated this new 

frequency 

A6 Option variable (See Section 3) 



OLDFCD 



16 Frequency and number of station in existing 

situation, to be changed to NFRE (NUMNEW 
cards) 
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Appendix III 
Subroutine MFCALC 



Transmitter cards 

Lat. Long. Freq. 

Prog. Source 



V 



Great circle 

distance between 

transmitters. 



Number of 'hops' 

and distance 

per 'hop' 



V 



Take-off angle 

and v.r.p. correction 

(in accordance with 

aerial type) 



V 



Field strength, 

corrected for 

v.r.p. and 

transmitter power 



V 



Frequency 
difference 



V 




Weighting curve 

protection ratio 

(different 

programmes) 
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